OPEN 3 ACCESS Freely available online 



•0-PLOS I o-^E 



Association between Polymorphisms in the TSHR Gene 
and Graves' Orbitopathy 

clKkhM updates 

Beata Jurecka-Lubieniecka^% Rafal Ploski^, Dorota Kula\ Konrad Szymanski^, Tomasz Bednarczuk^, 
Urszula Ambroziak^ Kornelia Hasse-Lazar\ Lidia Hyla-Klekot^, Andrzej Tukiendorf^, Zofia Kolosza^, 
Barbara Jarzab^ 

1 Maria Sklodowska-Curie Memorial Cancer Center and Institute of Oncology, Gliwice Branch, Department of Nuclear Medicine and Endocrine Oncology, Gliwice, Poland, 

2 Department of Medical Genetics, College of Forensic Medicine, Medical University of Warsaw, Warsaw, Poland, 3 Department of Internal Medicine and Endocrinology, 
Medical University of Warsaw, Warsaw, Poland, 4 Pediatric and Oncology Center, Chorzow, Poland, 5 Department of Epidemiology, Gliwice Branch, Maria Sklodowska- 
Curie Memorial Cancer Center and Institute of Oncology, Gliwice, Poland 



Abstract 

Background: Graves' orbitopathy (GO) as well as Graves' disease (GD) hyperthyroidism originate from an autoimmune 
reaction against the common auto-antigen, thyroid-stimulating hormone receptor (TSHR). GO phenotype is associated with 
environmental risk factors, mainly nicotinism, as well as genetic risk factors which initiate an immunologic reaction. In some 
patients GO is observed before diagnosis of GD hyperthyroidism, while it can also be observed far after diagnosis. The 
intensity of GO symptoms varies greatly in these patients. Thus, the pathogenesis of GD and GO may correlate with different 
genetic backgrounds, which has been confirmed by studies of correlations between GO and polymorphisms in cytokines 
involved in orbit inflammation. The aim of our analysis was to assess genetic predisposition to GO in young patients (age of 
diagnosis <30 years of age), for whom environmental effects had less time to influence outcomes than in adults. 

Methods: 768 GD patients were included in the study. 359 of them had clinically evident orbitopathy (NOSPECS >2). 
Patients were stratified by age at diagnosis. Association analyses were performed for genes with a known influence on 
development of GD - TSHR, HLA-DRBl, cytotoxic T-lymphocyte antigen 4 (CTLA4) and lymphoid protein tyrosine 
phosphatase (PTPN22). 

/feii//fs;The rsl 79247 TSHR polymorphism was associated with GO in young patients only. In young GO-free patients, allele 
A was statistically more frequent and homozygous carriers had a considerable lower risk of disease incidence than patients 
with AG or GG genotypes. Those differences were not found in either elderly patients or the group analyzed as a whole. 

Conclusions: IK\\e\e A of the rsl 79247 polymorphism in the TSHR gene is associated with lower risk of GO in young GD 
patients. 
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Introduction 

Orbitopathy is an inflammatory autoimmune reaction involving 
the orbit [1-2]. The primary autoantigen in GO and GD is TSHR 
which plays a key role in triggering the onset of disease [3]. 
Autoimmune reactions against TSHR located on fibroblasts 
prompt changes in the orbit connective tissue, which create 
fundamental clinical signs of disease [1-2], [4], Despite a common 
autoantigen, clinical symptoms of orbitopathy are diverse, starting 
from minimal symptoms to severe sight-threatening forms. 
Approximately one half of Graves' patients do not have clinically 
relevant GO [5]. GD hyperthyroidism may coexist with GO, 
however, most typically GO precedes or follows GD hyperthy- 
roidism at intervals ranging from a few months to a few years [6] . 



Taking into account the diverse nature of tire GO phenotype, 
one may presume the etiological factors behind GD are different 
from those leading to GO development. It has been proven that 
GO is a result of interactions between both genetic and 
environmental factors [7-9]. It has been confirmed that GO risk 
increases with age and is considerably higher in smokers [7], [10— 
11]. 

However, the role of genetic factors in the development of GO 
remains controversial. 

Current studies suggest a correlation between polymorphisms in 
cytokines involved in orbit inflammation and GO, but not GD 
[12-16]. However, a majority of studies do not find differences in 
polymorphisms within protein-coding genes involved in immune 
responses or the TSHR gene, between GD patients with and 
without orbitopathy [9], [17-20]. 
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Genetic factors may play a greater role than environmental 
factors in young GD patients [21-23]. An understanding of 
genetic predisposition to GO in young patients would allow the 
introduction of preventative measures, including smoking avoid- 
ance as well as early detection and treatment of hyperthyroidism 
[24] . In light of currently insufficient methods to treat GO, the 
potential for prevention is of particular importance [5] . Thus our 
study sought to analyse the genetic background of orbitopathy in 
different age groups, with special emphasis on genes with a known 
role in development of GD: TSHR, HLADRBl, CTLA4 and 
PTPN22. 

The study revealed an association between an TSHR polymor- 
phism and orbitopathy development in a young PoKsh population. 

Methods 

Genomic DNA was extracted from whole blood according to 
standard protocols by salt protein precipitation. HLA-DRBl 
typing was performed using sequence-specific oligonucleotides 
(SSO, Innolipa HLA-DRBl, Innogenetics, Gent, Belgium) and 
sequence-specific primers (MSSP Class II DRB Only, One 
Lambda, Dynal All Set SSP DR test, Dynal Biotech, Oslo, 
Norway) (Kula 2006) [25] in the Gliwice study, and by Dynal All 
Set SSP Dr test (Dynal Biotech, Bromborough, Wirral, UK) in the 
Warsaw study (Bednarczuk 2004) [26]. CTLA-4 and PTPN22 
were analysed by PCF-RFLP methods, as previously described 
Jurecka-Lubieniecka 2013 [27], Kula 2006 [25], Bednarczuk 
2003 [28], Skorka 2005 [29]) [Table 1]. In the Gliwice smdy PGR 
was performed with 0.5 units Hot Star Taq polymerase (Qiagen). 
PGR products were visualized on a 2% agarose gel stained with 
ethidium bromide and digested with the appropriate restriction 
endonuclease for 3 h at 37°C. Digestion products were separated 
on a 3% agarose gel. Genotyping of rsl79247 and rsl2101255 
polymorphisms in the TSHR gene was performed using TaqMan 
SNP genotyping technology (AppKed Biosystems, Foster City, 
USA) according to the manufacturer's protocol (Ploski 2010) [30]. 

Patients 

The analysis included 768 GD patients. Patients were divided 
into two groups using the criterion of age of onset: younger 
patients being diagnosed at ^30 years of age (n = 226), and older 
patients with an age of onset >30 years of age (n = 542). Patients 
were consecutively recruited in the Department of Nuclear 
Medicine and Endocrine Oncology, Centre of Oncology in 
Gliwice, Poland (n = 370) and in the Department of Endocrinol- 



ogy, Medical University of Warsaw, Poland (n = 398). AH GD 
patients in the cohorts investigated in the current study were 
unrelated and gave informed written consent. The project was 
approved by the relevant local research committees, including 
Medical Research Centre, Polish Academy of Science ethics 
committee and the Maria Sklodowska-Curie Memorial Cancer 
Centre and Institute of Oncology, Gliwice Branch ethics 
committee. Analysis was performed between 2003 and 2011, 
and median time of patient observation was approximately 5 
years. AH individuals were Caucasian. Criteria for the diagnosis of 
GD were the same for both patient groups, and were based on 
clinical and biochemical symptoms of hyperthyroidism, diffuse 
thyroid radioiodine uptake and/or detectable TSHR autoanti- 
bodies (Brahms). The severity of orbitopathy was assessed 
according to the NOSPECS classification. The clinical character- 
istics of GD patients are presented in Table 2. 

Statistical analysis 

Statistical analysis was performed using the statistical program 
STATA12.0. To ensure genotyping accuracy, Hardy-Weinberg 
equilibrium was calculated. The distribution of genotypes and 
alleles were compared between groups using the Chi-square test 
and Fisher's exact test [where appropriate]. AH tests were two- 
sided, and p<0.05 was considered significant. Odds ratios [ORs] 
were calculated with 95% confidence intervals (95% CI) to 
measure the effect of study factors on risk of GO. Association 
analyses of polymorphisms were performed for age of diagnosis 
30 years vs. >30 years), severity of orbitopathy based on 
NOSPECS classification (NOSPECS 0-1 vs. NOSPECS 2-6), 
and smoking status. A multivariate logistic regression was used to 
determine the independent association between GO and the above 
factors. 

To establish types of patients, a cluster analysis was applied. 
This analysis formed subsets of patients (clusters) by minimizing 
within-group variance and maximizing between-group variance 
using the chosen measure of dis/simUarity. The following risk 
factors were taken into account: smoking, TSHR rs 179247 
polymorphism, and age at GO onset. Once the patient groups 
were structured in clusters, they were able to be interpreted to 
explain the membership of patients in the clusters and their 
underlying predicting factors. This method may provide an 
alternative explorative platform for the identification of such 
predictors. Taxonomy uses a wide range of algorithms to 
determine the distance between objects. In our study, we used 



Table 1. Polymorphisms and methods used in the Gliwice study. 







Gene 


Polymorphism (formerly) 
rs (HapMap) 


Location 


Method 


PCR Primer F 5 -3' 
Primer R 5 -3' 


PCR Annealing 
temperature 


RFLP 
enzyme 


HLA-DRBl 




exon 2 


PCR-SSP/PCR-SSO 








CTLA-4 


A(49)G rs231775/rs57563726 


codon 1, 
exon 3 


PCR-RFLP 


F:CCAAGTCTCCA CTTAGTTATCC 
R:CCTCCATCTTC ATGCTCC 


55,1 °C 


Bst71 1, New 

England 

Biolabs 


PTPN22 


C(1858)T rs2476601 


codon 620, 
exon 14 


PCR-RFLP 


FiTCACCAGCTTC CTCAACCACA 
R:GATAATGTTGC TTCAACGGAATTT 


60°C 


Xcml, New 

England 

Biolabs 


TSHR 


rsl 79247 


intron 1 


TaqMan SNP 
genotyping 








TSHR 


rsl2101255 


intron 1 


TaqMan SNP 
genotyping 









doi:l 0.1 371 /journal.pone.Ol 02653.t001 
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Table 2. Clinical characteristics of patients with GD, N = 768. 





Gender (female: male) 


617:151 


Age at diagnosis of GD in years (mean ± SD) 


40.3 ±14.49 


GO present (N0SPECs=:2) 


359 (46.7%) 


Tobacco smokers 


322 (41.9%) 


Disease duration in years: (mean ± SD) 


2.72±4.38 


doi:l 0.1 371 /journal.pone.Ol 02653.t002 



the Marczewski-Steinhaus algorithm, which reUes on the use of a 
symmetric difference between objects [31]. 

Results 

Allele frequency and genotype distribution analysis 

The genotype frequency in the GD group was in Hardy- 
Weinberg equilibrium. 

In order to assess genetic predisposition to GO, we compared 
patterns of polymorphisms in the TSHR, PTPN, CTLA4 and 
HLA DRBl genes in patients with and without orbitopathy. 
Results were analysed both for the group as a whole group, as well 
as for subgroups of younger (age at diagnosis 30 or less) and older 
(age at diagnosis greater than 30) patients. 

When the group was analyzed as a whole, as well as when older 
patients were analyzed alone, there was no difference found in 
allele frequency or distribution of genotypes for any of the 
analyzed polymorphisms. Analysis of the younger patient group 
revealed significant differences in the presence of polymorphisms 
in relation to development of GO. Results of allele and genotype 
frequency analyses are presented in Figure 1. 

TSH receptor. Two polymorphisms in the TSHR gene 
(rs 179247 and rsl2101255) with a known association to GD in 
Polish populations were analyzed. It has been proven that allele A 
of rs 179247 and allele T of rs 12 10 1255 are more frequent in GD 
patients. Allele G of rsl79247 and allele C of rsl2101255 are more 
common in the healthy population [30]. In our study, in the group 
of younger patients, allele A of polymorphism rsl79247 was found 
significantly more frequently in patients without GO in compar- 
ison to those with GO (58.2% vs. 44.5%, p = 0,013, OR = 0.58) 
(Figure 1). In this group there was a statistically significant 
difference in genotype distribution (p = 0.038). The presence of a 
homozygous AA was associated with a significant reduction in risk 
of disease incidence, as compared to patients with AG or GG 
genotypes (p = 0.019, OR = 0.43) (Figure 1). Analysis of the 
rsl2101255 polymorphism in younger patients showed an 
increased presence of allele T and homozygous TT genotypes in 
patients without orbitopathy. This increase in frequency fell short 
of statistical significance (42.3% vs. 32.2% p = 0.059 and 20.2% vs. 
9.6%, p = 0.063 respectively). Of the 599 patients with GD, in 
which the TSHR polymorphisms were determined, analysis of 1 77 
patients young with and without GO showed a relationship 
between polymorphism rs 179247 of TSHR gene and the 
occurrence of GO (p = 0.013) with the power of 76%. 

CTLA4. Allele G was found more frequently in younger 
patients with GO in comparison to those without GO. However, 
this observation was not statistically significant (59.8% vs. 48.2%; 
p = 0.075; OR= 1.60) (Figure 1) 

PTPN22. Arrangement of alleles and genotypes of PTPN22 
C(1858)T/ rs2476601 did not differ in patients with and without 
GO (Figure 1). 



HLA DRB1*03. Arrangement of alleles and genotypes of 
polymorphism HLA DRB1*03 did not differ in patients with and 
without GO (Figure 1). 

Haplotype frequency analysis 

The analyzed polymorphisms, rsl2101255 and rsl79247 show 
linkage disequilibrium (LD). D'is very close to 1.0 (D' = 0.97), 
which indicates a strong linkage disequilibrium at p = 2.2204e-16. 
The most common haplotype in patients with GD is the haplotype 
comprising two AT-risk alleles. It occurs more frequently in 
patients with GD than in healthy subjects (37.6% vs. 26.2%, OR 
= 1.72, p = 0.00002) (data unpublished). 

Logistic regression 

Multivariate logistic regression was performed using GO as the 
dependent variable and the studied polymorphisms, age of onset 
and smoking status as the independent variables. 

Results of the logistic regression confirmed observations from 
analyses of alleles and genotypes. Genetic predisposition and 
smoking are independent risk factors with influences on develop- 
ment of orbitopathy in younger patients. The presence of the 
rsl 79247 polymorphism in the TSH receptor gene significantly 
lowers the risk of GO incidence (OR = 0.571, p = 0.012). 
However, in older patients smoking is the only significant factor 
influencing development of orbitopadiy (OR= 1.802, p = 0.001) 
(Table 3). 

Taxonomy 

In order to confirm our results, additional analyses of GO 
incidence among GD patients were performed and compared with 
results of the GO risk factor analysis (rsl 79247 polymorphism, 
smoking and age at diagnosis) obtained using logistic regression. 
Taxonomy was used to group patients in terms of similar values for 
chosen parameters. The group was divided into three subgroups 
(Figure 2): 

Subgroup A included older smoking patients (median age 43), 
among which all rsl 79247 genotypes (AA, AG, GG) were found. 
Polymorphism frequency was not associated with the presence or 
absence of GO. 

Subgroup B was comprised of non-smoking patients with a GG 
genotype. 

The most interesting to us was subgroup C, which included 
allele A carriers, non-smokers and had the lowest median age (35 
years). In this subgroup a significant discrepancy in genotype 
distribution (AG and AA) was revealed depending on age of 
diagnosis and orbitopathy occurrence. AUele A carriers were 
dominant in the group of patients without orbitopathy (N = 145 vs. 
N = 90) (data unpublished). The difference was more significant in 
young (age of diagnosis ^30 years of age) allele A carriers (N = 64 
vs. N = 28) (data unpublished). Whereas orbitopathy was not 
present in 86% of young AA homozygote carriers (subgroup C3) 
(N = 30 vs. N = 5) (Figure 3). 
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gen/SNP 


allele/ 
genotype 


all patients 


patients < 30 years at diagnosis 


GD 
without 
GO 


GDwith 
GO 


p value 

(OR 
95%CI) 


GD 
without 
GO 


GD w'lXYi 
GO 


p value 
(OR) 


TSHR 
rsl79247 




N-316 
(100%) 


N-283 
(100%) 




N- 104 
(100%) 


N= 73 
(100%) 






302 
(47.8%) 


280 
(49.5%) 


0.56 
(0.93, 0.74 
1.18) 


87 
(41.8%) 


81 
(55.5%) 


0.013 
(0.58, 0.37 
0.90) 


A 


330 


286 


121 


65 


GG 


77 
(24 4%) 


77 


0.71 


20 

CIO ''%^ 


22 
(30 1%) 


0.038 


AG 


148 
(46.8%) 


126 

(44.5%) 


47 
(45.2%) 


37 
(50.7%) 


AA 


(28.8%) 


(28.3%) 


(35.6%) 


(19.2%) 


AA vs 
GG+AG 






0.88 
(0.97, 0.68 
1.39) 






0.019 
(0.43, 0.20 
0.91) 


TSHR 
rsl2101255 




N = 314 
(100%) 


N = 282 
(100%) 




N= 104 
(100%) 


M = 73 
(100%) 






383 
(61.0%) 


352 
(62.4%) 


0.81 
(0.96, 0,75 
1.24) 


120 
(57.7%) 


99 
(67.8%) 


0.059 
(0.65, 0.40 
1.03) 


T 


245 
f n%i 


212 


88 


47 


CC 


119 
(37.9%) 


112 

(39.7%) 


0.96 


37 
(35.6%) 


33 
(45.2%) 


0.131 


CT 


(46.2%) 


(45.4%) 


(44.2%) 


(45.2%) 


TT 


50 
(15.9%) 


42 
(14.9%) 


21 
(20.2%) 


7 (9.6 %) 


TT vs 
CT+CC 






0.93 
(0.95,0.58 
1 55'1 






0.063 
(0.42, 0.14 
1 1 n 


CTLA4 
A(49)G 
rs231775 




N . 325 
(100%) 


N . 304 
(100%) 




N .83 
(100%) 


N =56 
(100%) 




A 


335 
(51.5%) 


297 
(48.8%) 


0.37 
(1.11, 0.89 
1.40) 


86 (51.8%) 


45 (40.2%) 


0.075 
(1.60, 0.96 
2.68) 


G 


315 
(48.5%) 


311 
(51.2%) 


80 (48.2%) 


67 (59.8%) 


AA 


88 (27.1%) 


75 (24.7%) 


0.64 


25 (30.1%) 


10 (17.9%) 


0.19 


AG 


159 
(48.9%) 


147 
(48.4%) 


36 (43.4%) 


25 (44.6%) 


GG 


7R OA fWA 


ft7 ill n%\ 


22 (26 5%) 


")! n7 '^%\ 


PTPN22 
C(1858)T 
rs2476601 




N = 298 
(100%) 


N = 261 
(100%) 




N = 90 
(100%) 


N = 52 
(100%) 




C 


486 
(81.5%) 


430 
(82.4%) 


0.78 
(0.95, 0.69 
1.30) 


143 
(79.4%) 


77 (74.0) 


0.37 
(1.35, 0.73 
2.48) 


T 


110 

(18.5%) 


92 (17.6%) 


37 (20.6%) 


27 (26.0%) 


CC 


194 
(65.1%) 


180 
(69.0%) 


0.12 


57 (63.3%) 


31 (59.6%) 


0.30 


CT 


98 (32.9%) 


70 (26.8%) 


29 (32.2%) 


15 (28.8%) 


TT 


6 (2.0%) 


11 (4.2%) 


4 (4.4%) 


6 (11.5%) 


HLA-DRBl 
HLADRB1»03 




N . 209 
(100%) 


N = 220 
(100%) 




N = 58 
(100%) 


N =42 
(100%) 




DR3- 


317 
(75.8%) 


343 
(78.0%) 


0.51 
(0.89, 0.64 
1.24) 


83 (71.6%) 


55 (65.5%) 


0.45 

(1.33, 0.69 
2.53) 


DR3 + 


101 
(24.2%) 


97 (22.0%) 


33 (28.4%) 


29 (34.5%) 


DR3 -/- 


113 
(54.1%) 


128 
(58.2%) 


0.69 


26 (44.8%) 


15 (35.7%) 


0.50 


DR3 -/+ 


91 (43.5%) 


87 (39.6%) 


31 (53.4%) 


25 (59.5%) 


DR3 +/+ 


5 (2,4%) 


5 (2.3%) 


1 (1.7%) 


2 (4.8%) 



Figure 1. Frequencies of alleles and genotypes in Graves' patients with and without orbitopathy. 

doi:10.1371/journal.pone.0102653.g001 
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Table 3. Results of multiple linear regression analysis. 







OR 


95% confidence interval 


p-value 


Age at diagnosis £30 TSHR rsl 79247 Smoking 


0.571 


0.37-0.88 


0.012 




1.879 


0.97-3.66 


0.063 


Age at diagnosis >30 Smoking 


1.802 


1 .27-2.56 


0.001 



doi:l 0.1 371/journal.pone.0102653.t003 



Discussion 

The TSHR gene encodes the thyroid-stimulating hormone 
receptor (TSHR) - a protein that plays a role in the pathomech- 
anism of GO. Thus, polymorphisms in this gene may result in a 
predisposition for GO [1-3], [9]. TSH receptor is the dominant 
autoantigen involved in the development of an immune response 
carried out by immunocompetent cells expressing HLA class II 
complex on the surface [32] .While the cascade of immunological 
response is generally well known, the mechanisms leading to the 
origin of the MHC II-TSHR complex need to be analyzed more 
extensively [1], [4], [9]. The etiology of disease may result from 
genetic differences in the TSHR gene, leading to TSHR protein 
structural dififerences [33-34]. Intronic polymorphisms - not 
exonic - are responsible for generation of different receptors forms 
[34]. 

The TSHR gene association with GD was well established. Two 
independent studies identified strong association within TSHR 
intron 7 in Japanese and TSHR intron 1 in UK Caucasian 
cohorts, which provided the first convincing evidence for TSHR 
association with GD [35-36]. Recent detailed association mapping 



of 98 SNPs across 800 Kb of the TSHR region, refined association 
to within 40 Kb of TSHR intron 1 in 768 GD patients and 768 
controls of UK Caucasian origin. Strongest SNP associations were 
rsl79247 (OR= 1.53) and rsl2101255 (OR= 1.55), which were 
replicated in Polish-Caucasian cohorts. A different study showed 
that two SNPs (rsl 2 10 1261, rsl 79243) in intron 1 are associated 
with GD in Chinese population [37]. However these polymor- 
phisms are in linkage disequilibrium with SNP's previously 
reported in Caucasian populations. These results strongly suggest 
that SNP's in intron 1 of TSHR are risk factors for GD. The role 
of SNP's in intron 7 remain to be established. 

In our work we analyzed the association of SNPs rsl2101255 
and rsl 79247 with orbitopathy as their relationship with the 
development of GD was proved in the Pohsh population [30] . 

The location of rsl 79247 and rsl2101255 within the first 
intron, close to the promoter region and start codon, may 
influence gene expression or posttranslational processes. Structur- 
ally and functionally altered TSHR isoforms become an 
immunotarget in patients with a genetic predisposition to GD. 
Polymorphisms rsl79247 and rsl2101255 were associated with 




Figure 2. Dendrogram: age at diagnosis, smoking, TSHR rsl 79247. Group A: N = 249, median age 43 years, smokers, genotypes: GG, AG, AA, 
with GO: N = 138 (55%), without GO: N = 112 (45%). Group B: N = 73, median age 41 years, non smokers, genotypes: GG with GO: N = 33 (45%), 
without GO: N = 41 (55%). Group C: N = 235, median age 35 years, non smokers, genotypes: AA,AG. with GO: N = 90 (38%), without GO: N = 145 (62%). 
Subgroup CI: N = 85, median age: 48 years, genotypes: AG with GO: N = 37 (44%), without GO: N=48 (56%). Subgroup C2: N = 59, median age: 22 
years, genotypes: AG With GO: N = 23 (39%), without GO: N = 36 (61%). Subgroup C3: N = 91, median age 35 years, genotypes: AA With GO: N = 30 
(33%), without GO: N = 61 (67%). 
doi:1 0.1 371/journal.pone.01 02653.g002 
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Figure 3. Incidence of GO in subgroup C3 with AA genotype TSHR rs 179247. 

doi:1 0.1 371 /journal.pone.01 02653.g003 



reduced expression of full length TSHR mRNA relative to two 
truncated splice variants in thyroid tissue [2], [4], [33-34]. 

Our analysis demonstrated an association between orbitopathy 
and the TSHR gene rs 179247 polymorphism. In younger patients 
without orbitopathy, in whom GD was diagnosed before 30 years 
of age, the A allele occurred significandy more often. Presence of 
an AA homozygous locus was associated with a significant 
reduction in risk of GO incidence. Independent taxonomy also 
confirmed the relation between TSHR polymorphism rs 179247 
and GO in younger patients. Orbitopathy was not present in 86% 
of young AA carriers (N= 30 vs. N = 5) (Figure 3). Patients with 
allele T or homozygous TT at the rsl2101255 polymorphism in 
TSHR demonstrated an increased frequency in GO incidence 
within this group of patients (Figure 1). The AT haplotype occurs 
most often in patients with GD, and is associated with an increased 
risk of disease (OR= 1.72, p = 0.0002). 

However, when all patients were analysed, or when only older 
patients (age at diagnosis >30) were analysed, the frequency of 
alleles present and genotype distributions of both polymorphisms 
did not differ in patients with or without GO. 

These findings showed an association between the presence of 
TSHR polymorphisms in younger patients and the incidence of 
orbitopathy. Allele A of rs 179247 was more frequent in GD 
patients regardless of age than in healthy persons [30]. In our 
study this allele showed a correlation with lesser risk of GO in 
young GD patients. 

This was also confirmed by preliminaiy epidemiologic studies 
showing that GO incidence appears in two peaks throughout life - 
after 40 years of age and after 60 years. In these patients 
environmental factors (mainly nicotinism) were the factors most 
closely correlated with illness occurrence [10], [24]. In younger 
patients GD-associated eye symptoms were mild and character- 
ized by a self-limiting course without involvement of the 
oculomotor nerve [38-41]. Current pubhcations have not revealed 
a correlation between TSHR polymorphisms and orbitopathy 
[19-20]. These .studies have focused on an rs2268458 polymor- 
phism localized within the first mtron of the TSHR gene. The 
development of GO showed no association with the TSHR 
polymorphism, however the patient group was relatively small and 
the median age was 48 [19]. In a second study, extended to 199 
patients, the population was dominated by patients with severe 
orbitopathy and orbital decompression, with a median age of 52. 
There were no differences revealed in the frequency of allele C 
amongst GD patients with or without GO [20]. Our analysis of 
older patients (greater than 30 years old at diagnosis) as well as the 



group analyzed as a whole (median age 40) yielded similar results. 
Taxonomic classification allowed the separation of a group of 
older smoking patients (group A, median age 45), in which there 
was no association between the rs 179247 TSHR polymorphism 
and orbitopathy (Figure 2). 

Syed et al. analyzed the influence of the rs2268458 TSHR 
polymorphism, as well as seven PTPN12 polymorphisms, on the 
phenotypes of GD. Association of the TSHR rs2268458 
polymorphism with mild/moderate orbitopathy was confirmed, 
however in interaction with PTPN12 polymorphisms only. 
Authors showed independent correlations of three PTPN12 
polymorphisms with orbitopathy and interaction of each with 
the rs2268458 TSHR polymorphism. No age related dependen- 
cies were found [42]. 

The role of other polymorphisms in development of orbitopathy 
remains unclear. 

The first studies looking at genetic background in the context of 
orbitopathy indicated an association between the A(49)G poly- 
morphism in the CTLA4 gene and an increased risk of GO [43]. 
These conclusions were confirmed in a later analysis of the same 
group [44] . Other studies of Caucasian groups did not confirm the 
relationship between the A(49)G polymorphism and GO [10], 
[45-47]. 

The HLADRB1*03 polymorphism occurs most commonly in 
GD patients of Caucasian origin. Some papers have confirmed an 
association between HLADRB1*03 and GO [48-49]. However, 
other authors have not found differences in the distribution of the 
HLADRB1*03 allele in GD patients with or without GO [10], 
[27], [50-51]. There is no sufficient evidence to support a 
correlation between the C(1858)T polymorphism in the PTPN22 
gene and GO [29], [42]. 

Our analysis of A(49)G polymorphisms in the CTLA4 gene and 
C(1858)T polymorphisms in the PTPN and HLADRB1*03 genes 
did not reveal any statistically significant discrepancies in allele 
frequencies or distribution in patients with or without orbitopathy 
in aU age groups. However, there was an observed tendency 
showing a more frequent occurrence of allele T of the C(1858)T 
polymorphism in the PTPN gene, the HLADRB1*03 allele and 
the G allele of the A(49)G polymorphism in the CTLA4 gene in 
young patients with GO. One limitation of this study was the 
inclusion of a relatively small group of patients in particular 
subgroups, which causes statistical difficulties and may indicate 
caution in the interpretation of negative results. Future studies wiU 
include a larger group of patients. 
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Difficulties in die assessment of genetic predisposition to GO are 
in large part due to genotype heterogeneity in autoimmunological 
disorders. Interactions between many genes with modest increases 
in relative risk influence the development of particular GD 
phenotypes, including orbitopathy [52]. Thus it is very important 
to precisely classify GO patients according to the clinical course of 
disease [18]. Analysis of genetic predisposition is also hampered by 
the necessity for long-term observation, due to undefined time 
frames for development of hyperthyroidism and GO, as well as the 
confounding effects of environmental factors [6], [53]. 

Introduction of the described age-grouping method may allow a 
simplification in assessment of genetic predisposition to GO 
because age grouping decreases the burden of long-term 
observation. In the young group of patients environmental 
influences have had less time to play a role in disease course, 
therefore these environmental effects are relatively minimized. 
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Lesser influence of the environment is an important factor to be 
considered in favor of the validity of a method dividing analyzed 
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To continue our study we have planned a wider assessment of 
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In summary, we have found that allele A of the rs 179247 
polymorphism in the TSHR gene is associated with lower risk of 
GO in young patients with GD. 

Author Contributions 

Conceived and designed the experiments: BJL AT BJ. Performed the 
experiments: BJL DK RP TB KS UA. Analyzed the data: AT ZK. 
Contributed reagents/materials/analysis tools: AT KHL LHK. Wrote the 
paper: BJL DK AT. 



21. Tomer Y, Menconi F, Davies TF, Barbesino G, Rocchi R, et al (2007) Dissecting 
genetic heterogeneity in autoimmune thyroid diseases by subset analysis. 
J Autoimmun 29: 69-77. 

22. Hemminki K, Li X, Sundquist J, Sundquist K (2010) The epidemiology of 
Graves' disease: evidence of a genetic and an environmental contribution. 
J Autoimmun 34: J307-13. 

23. Segni M, WoodJ, Pucarelli I, Toscano V, Toseano R, et al (2001) Clustering of 
autoimmune thyroid diseases in children and adolescent: a study of 66 families. 
J Pediatr Endocrinol 14 Supl 5: 1271-5; discussion 1297-8. 

24. Wiersinga WM, Bartalena L (2002) Epidemiology and prevention of Graves' 
orbitopatiiy. Thyroid 12(10): 855-60. 

25. Kula D, Bednarczuk T, Jurecka-Lubieniecka B, Polanska J, Jarzab M, et al 
(2006) Interaction of HLA-DRBl alleles with CTLA-4 in the predisposition to 
Graves' disease: the impact of DRB1*07. Thyroid 16(5): 447-53. 

26. Bednarczuk T, Hiromatsu Y, Seki N, Ploski R, Fukutani T, et al (2004) 
Association of tumor necrosis factor and human leukocyte antigen DRB 1 alleles 
with Graves' orbitopathy. Hum Immunol 65(6):632— 9. 

27. Jurecka-Lubieniecka B, Ploski R, Kula D, Krol A, Bednarczuk T, et al (2013) 
Association between age at diagnosis of Ciraves' disease and variants in genes 
involved in immune response. PLoS One 8(3):e59349. doi: 10.1371/journal. 
pone.0059349. Epub 2013 Mar 27. 

28. Bednarczuk T, Hiromatsu Y, Fukutairi T, Jazdzewski K, Miskiewicz P, et al 
(2003) Association of cytotoxic T-lymphocyte-associated antigen-4 (CTLA-4) 
gene polymorphism and non-genetic factors with Graves' orbitopathy in 
European andjapanese populations. EurJ Endocrinol 148(1):13— 8. 

29. Skorka A, Bednarczuk T, Bar-Andziak E, Nauman J, Ploski R, et al (2005) 
Lymphoid tyrosine phosphatase (PTPN22/LYP) variant and Graves' disease in 
a Polish population: association and gene dose-dependent correlation with age of 
onset. Clin Endocrinol (Oxf) 62(6):679-82. 

30. Ploski R, Brand OJ, Jurecka-Lubieniecka B, Franaszczyk M, Kula D, et al (20 1 0) 
Thyroid stimulating hormone receptor (TSHR) Intron 1 variants are major risk 
factors for Graves' disease inthree european Caucasian cohorts. PLoS One 
25;5(ll):el5512. doi: 10.1371/joumal.pone.0015512 

31. Tukiendorf A, Kazmierski R, Michalak S (2013) The taxonomy statistic 
uncovers novel clinical patterns in a population of ischemic stroke patients. PLoS 
ONE 8(7): e69816. doi:10.1371/journal.pone.0069816 

32. Sawai Y, DcCiroot L (2000) Binding ol human thyrotropin receptor peptides to a 
Graves' disease predisponsiiig human leukocyte antygen ealss 11 molecule. J Clin 
Endocrinol Mctab 85: 1176-9. 

33. Brand OJ, Cough SGL (2010) Genetics and thyroid autoimmunity and the role 
of tiie TSHR. Mol Cell Endocrinol 322: 135^3. 

34. Davies TF, Yin X, Latif R (2010) The genetics of the thyroid stimulating 
hormone receptor: history and relevance. Thyroid 20(7): 727—35. 

35. Hiratani H, Bowden DW, Ikegami S, Shirasawa S, Shimizu A, et al. (2005) 
Multiple SNPs in intron 7 of thyrotropin receptor arc associated with Graves' 
disease. J Clin Endocrinol Metab 90: 2898-2903. 

36. Brand OJ, Barrett JC, Simmonds MJ, xNewby PR CJ, McCabe CJ, et al (2009) 
Association of the thyroid stimulating hormone receptor gene (TSHR) with 
Graves' disease. Hum Alolceular Gen 18: 1704 1713. 

37. Liu BL, Yang SY, Liu W, Xue LC^ Chen -X, el al :201 1; Refined assotiation of 
TSH receptor susceptibility locus to Graves' disease in the Chines Han 
population. EurJ Endocrinol 170: 109-119. 

38. Goldstein SM, Katowitz WR, Moshang T, KatowitzJA (2008) Pediatric tiiyroid- 
associated orbitopathy: the children's hospital of Philadelphia experience and 
literatiire review. Thyroid 18(9): 997-9. 

39. Chan W, Wong GWK, Fan DSP, Cheng AC, Lam DS, et al (2002) 
Ophthalmopathy in childhood Graves' disease. Br J Ophthalmol 86(7): 740—2. 

40. Krassas GE (2004) Treatment of juvenile Graves' disease and its ophthalmic 
complication: the 'European way'. Europ J Endocrinol 150: 407-14. 



PLOS ONE I www.plosone.org 



7 



July 2014 I Volume 9 | Issue 7 | e102653 



The RTSH Gene and Graves' Orbitopathy 



41. Krassas (iE, Scgni M, Wicrsinga WM (2005) Childhood Ci-ravcs' orbitopathy: 
results of a European questionnaire study. Eur J Endoerinol 153: 515—20. 

42. Syed AA, Simmonds MJ, Brand OJ, Frankiyn JA, Gough SCL, et al (2007) 
Preliminary evidence for interaction of PTPN12 polymorphism with TSHR 
genotype and association with Graves' orbitopathy. Clin Endocrinol 67: 663—7. 

43. Vaidya B, Imrie H, Perros P, Dickinson J, McCarthy MI, et al (1999) Cytotoxic 
T lymphocyte antygen-4 (CTLA-4) gene polymorphism confers susceptibility to 
thyroid associated orbitipathy. Lancet 354: 743—4. 

44. Vaidya B, Oakes EJ, Imrie H, Dickinson AJ, Perros P, et al (2003) C.rrLA4 gene 
and Graves' disease: association of Graves' disease with the CrrLA4 exon 1 and 
intron 1 polymorphisms, but not with the promoter polymorphism. Clin 
Endocrinol (Oxf) 58: 732-5. 

45. AUahabadia A, Heward JM, Nithiyananthan R, Gibson SM, Reuser TT, et al 
(2001) MHC class 11 region, GTLi\4 gene, and orbitopathv in patients with 
Graves' disease. Lancet 358: 984-5. 

46. Bednarczuk T, Hiromatsu Y, Eukutani 'L, Jazdzewski K, Miskiewiez P, et al 
(2003)Assoeiation of cytotoxic 'L-lymphoeyte-associated antigen-4 ((.rrLi'\-4) 
gene polymorphism and non-genetic factors with Graves' orbitopathy in 
Europeaul andjapamese populations. Eur J Endocrinol 148: 13—8. 



47. Buzzetti R, Nistico L, Signore A, Caseino 1 (1999) CTLA-4 and HLA gene 
susceptibility to thyroid-associated orbitopathy. Lancet 354: 1824. 

48. Semana G, Mehra NK, LepetitJ, Horse J (1997) HLA and Graves' disease: 12th 
International Histocompatibility Workshop study. In Genetic Diversity of HLA 
Functional and Medical Implication, Medicad and Scientific International 
Pubhsher, EDC, Paris, France, 1, 407. 

49. Farid NR, Stone E, Johnson G (1980) Graves' disease and HLA, clinical and 
epidemiologic associations. CHn Endocrinol 13; 535^4. 

50. Kendall-Taylor P, Stephenson A, Stratton A, Papiha SS, Perros P, et al (1988) 
Differentiation of autoimmune orbitopathy from Graves' hyperthyroidism by 
analysis of genetic markers. Clin Endocrinol 28(6): 601-10. 

51. Wectman AP, So AK, Warner CA, Foroni L, Fells P, et a! (1988) 
Immunogenetics of (jraves' orbitopathy. Clin Endoerinol 28(6): 619—28. 

52- Davies TE (2007) Really significant genes for autoimmune thyroid disease do not 
exist - so how can wc predict disease. Thyroid 17: 1027-29. 

53. Brix TH, Kyvik KO, Christensen K, Hegediis L (2001) Evidence for a major 
role of heredity in Graves' disease: a population-based study of two Dguiish twin 
cohorts. J Clin Endocrinol Metab 86(2): 930-^. 



PLOS ONE I www.plosone.org 



8 



July 2014 I Volume 9 | Issue 7 | e102653 



